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Abstract 

Predictions for the spin dependent structure function g\ of the nucleon are presented. We use an unified 
approach incorporating the LO DGLAP evolution and the resummation of double logarithmic terms ln 2 (x). 
We show, that the singular input parametrisation as x — > can be a substitute of the ln 2 (x) resummation. 
An impact of the 'more running' coupling is discussed. We determine the contribution to the Bjorken sum 
rule solving the evolution equation for the truncated moment of gi S . A comparison with the re-analysed 
HERMES and COMPASS data is given. 
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1 Introduction 
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Experimental data confirm (at least for Q 2 > 1 GeV 2 ) the theoretical predictions of an increase of the 
nucleon structure functions at small values of the Bjorken x. The low- x behaviour of both spin averaged 
and spin dependent structure functions is controlled by the double logarithmic terms (a s ln 2 (x)) n In 
a unpolarised case, this singular PQCD behaviour is however overridden by the leading Regge contribution 
present in the input parametrisation gj. The situation is quite different in the spin-dependent case, where the 
double logarithmic effects are very important. The resummation of the ln 2 (x) terms at low x goes beyond the 
standard LO and NLO PQCD evolution of the parton densities. Double logarithmic contributions become 
VJ^ ' essential for x ~ 0.01, where there is little experimental data. Determination of the sum rules and the 

nucleon spin decomposition among partons requires knowledge of the structure functions over the entire 
\Q • region of the variable x £ (0; 1). Therefore the small- :r behaviour of the spin dependent parton distributions 

is a topic of the intensive theoretical investigations. Standard approach describing structure functions is 
based on the DGLAP-Q 2 evolution equation via two-step convolution: of the initial parton densities and 
splitting functions and then of the evolved parton distributions and the coefficient functions. Because there 
is no way to calculate the initial parton densities, which have a nonperturbative origin, they must be put 
'by hand'. Different parametrisations of the initial gluon and quark densities known in literature e.g. [3[B] 
are singular when x — > 0. This choice enables to study DIS phenomena within DGLAP approach not only 
for the large- x region but for the small- a; one as well. Singular terms ~ x~ x can be a substitute of the 
double logarithmic ln 2 (x) resummation, which is absent in the standard DGLAP scenario. This problem 
has been widely discussed and argued in [7]. Also important is the problem of the a 3 dependence of the 
QCD evolution. Following [8] we take into account the running coupling effects not via a s (Q 2 ) but with use 
of the more running a s (Q 2 jx). This approach is better justified at small values of x, whereas for large x ~ 1 
leads to the usual DGLAP coupling a s (Q 2 ). It seems to be reasonable to study an impact of the double 
logarithmic and running coupling effects on theoretical predictions for spin structure functions. 

In this paper we present the unified approach, in which the familiar Q 2 evolution is extended by the In 2 (x) 
resummation. In our analysis we use so-called unintegrated parton distributions and solve the combined 
LO DGLAP +ln 2 (x) evolution equation with help of the Chebyshev polynomial technique. We take into 
account the 'very running' coupling effects at small x and discuss the role, they play. We also show that the 
singular input parametrisation of the parton distributions can be some kind of substitution for the double 
logarithmic terms, missing in the standard DGLAP approximation. Our theoretical predictions for the spin 
dependent structure function g\ are compared with recently re-analysed HERMES and COMPASS data. 

The content of this paper is as follows. In Section 2 we recall the unified approach incorporating DGLAP 
evolution of structure functions and the double logarithmic a"ln 2n (x) terms, which are essential in the 
small- x region. Section 3 is devoted to the impact of running coupling effects on the gi results in the small- x 
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region. In Section 4 we show that the singular initial parton densities ~ x~ can mimic the resummation 
of double logarithmic terms ln 2 (x). Using this fact, in Section 5 we solve the evolution equation for the 
truncated moments themselves and obtain contribution to the Bjorken sum rule. We also present numerical 
predictions for the structure function gi and compare them to re-analysed HERMES and COMPASS data. 
Finally, in Section 6 we summarise our results. 



2 Unified ln 2 (x)+LO DGLAP approach 

The structure functions of the nucleon can be expressed in terms of the parton distributions. These depend 
on two kinematic variables: the Bjorken x and Q 2 = — q 2 with q being the four-momentum transfer in 
the deep-inelastic lepton-nucleon scattering (DIS). The scaling variable is defined as x = Q 2 /(2pq), where 
p is the nucleon four-momentum. The strong interactions between quarks and gluons cause the changes 
in the parton densities. For medium and large x, the evolution with Q 2 of the parton distributions is 
well described by the standard Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations [5]-|12|. This 
approach which effectively sums up the leading ln(Q 2 ) terms is however incomplete at small x, where another 
large logarithm - ln(l/x) becomes essential and which leading powers a" ln 2n (x) needs to be resummed. The 
double logarithmic terms ln 2 (x) come from the ladder diagrams with quark and gluon exchanges along the 
chain. Treating both potentially large logarithms ln(Q 2 ) and ln(l/x) on equal footing, the authors of |14]-[16] 
obtained equations which incorporate DGLAP evolution and ln 2 (x) terms as well. The double logarithmic 
effects go beyond the standard LO and even NLO Q 2 evolution of the spin dependent parton distributions 
and significantly modify the Regge pole model expectations for the structure functions. Theoretical analyses 
of the small- a; behaviour of the polarised structure functions |17| predict that resummation of the double 
logarithmic terms (a s ln 2 (x)) n leads to the singular form as x — > 0: 

9i S ' S (x,Q 2 ) ~ x~ Ajvs ' s , (1) 

where Ajvs ~ 0.4, As » 0.8 and gi S '" S denotes nonsinglet or singlet part of the polarised structure function 
of the proton. For larger but still low x € (10~ 5 ; 10~ 2 ), gi is less steep with the slope A ~ 0.2 — 0.3 for the 
nonsinglet part [TJ. This power-like behaviour x~ x remains significantly steeper than the DGLAP solution 
in absence of the singular input parametrisations of parton densities 



g? aLAP (x - 0) ~ exp lnln(QVA 2 QOT ) . (2) 

The unified equation, which includes the LO DGLAP evolution and ln 2 (x) terms resummation reads: 



f(x,Q 2 ) = f ( x ) + J^ g- AP®/(x,fc' ; 
Ql 



4 
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DGLAP 

l Q 2 /z 

f dz f dk a s tx , 2 \ 
J z J 2tt 1 U' fc ) 



LN 2 (X) LADDER 

+ Bremsstrahlung corrections , (3) 



LN 2 (X) NON LADDER 

where ® abbreviates a Mellin convolution over x 

i 



(AP ® /) (x, Q 2 ) = I — AP (- ) f(y, Q 2 ) , (4) 



y \y 

X 

AP denote the polarised version of the splitting function P and / is the unintegrated distribution, related 
to the ordinary polarised parton distribution Ap{x, Q 2 ) via 

2 _ 8Ap(x,Q 2 ) 

The double logarithmic terms come from ladder-type graphs as well as from the nonladder ones which 
represent radiative corrections [2], |18]-[20], In a case of the nonsiglet polarised structure functions the 
contribution of nonladder diagrams is negligible. However for the singlet spin dependent structure functions, 
besides the ladder graphs, one has to include Bremsstrahlung corrections [3], which are important. The 
full evolution equations for nonsinglet and singlet unintegrated parton distributions within DGLAP+ln 2 x 
approach have been presented in |14l 1 15] . This forms the basis of our analysis in the next section, where we 
discuss the modified running coupling effects at small x. 
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3 Running coupling a s in the small- £ region 



DGLAP formalism uses the following prescription for the running coupling (in the lowest order): 

a s = a s (Q ) = 



(33 — 2JV/) ln-j$?- 



(6) 



where Nf is the number of active quark flavours and Aqcd ~ 200 MeV is the QCD cut-off parameter. It 
has, however, been argued that in the small- a; region eq.© should be rearranged into the following form [5]: 



a s = a s (Q 2 /z 



(7) 



with z being the longitudinal momentum fraction of a parent parton, carried by a next generation parton. 
In this way a s becomes 'very running' i.e. runs in each ladder rung depending on the gluon virtuality. This 
prescription of a s , widely discussed also in 21, 22], has been used e.g. in [T] [23] within double logarithmic 
effect ln 2 (x) resummation. Here, we study an impact of the different a s parametrisation on the polarised 
parton densities. Because 

a g (Q 2 /z) < a s (Q 2 ) , (8) 

in the case of the 'very running' a s {7]), the growth of the parton distributions in the low- a; region is damped. 
A scale of the damping for nonsinglet and singlet (gluons) distributions is shown in Figs. 1-2, where we plot 
the ratio 

A P (a s (Q 2 )) 

Ap(a.(Q7*)) 1 ' 

as a function of x. Here, Ap denotes the nonsinglet (valence) Aq NS and the gluon AG distribution function 
respectively. One can see, that the difference becomes essential at x ~ 0.01 and the impact of the running 
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for the polarised nonsinglet quark distribution Aq NS = Au 
10 GeV 2 . Solid: unified DGLAP +ln 2 (x) approach, dotted: DGLAP alone. 



Figure 1: The ratio 



Ad as a function of x at 



coupling effects for the singlet case is larger than for the nonsinglet one. Double logarithmic resummation 
additionally amplifies the split between results in comparison to the pure DGLAP approach. At very 
small x = 10 -5 we find the ratio © about 2 for the nonsinglet polarised distribution and above 6 for the 
polarised gluons. Our estimations of R @ show that for the small values of Bjorken parameter x < 10~ 2 
the coupling a s (Q 2 ) should be replaced by a B (Q 2 /z). In standard DGLAP analysis, where rather large-a; 
region is considered, this modification converts into a s (Q 2 ) {z ~ 1). Parametrisation of the coupling a s is 
not the only crucial point in the low-x analysis of structure functions. Another problem are initial parton 
distributions at low Q 2 , ~ 1 GeV 2 , which enter into the evolution equations. The behaviour of the quark and 
gluon distributions at very small x is mainly generated by the double logarithmic ln 2 (x) effects. Therefore 
singular as x — > inputs ~ x~ ai seem to be needless in PQCD analysis, unless one does not consider 
ln 2 (x) terms. Within standard DGLAP approach, parametrisations in a form ~ x~ ai can be regarded as a 
substitute of the missing double logarithmic effects resummation. In the next section we discuss this problem 
in detail. 
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Figure 2: The ratio §§§ for the polarised gluon distribution AG as a function of x at Q 2 — 10 GeV 2 . Solid: 
unified DGLAP+/n 2 (x) approach, dotted: DGLAP alone. 

4 Singular input parametrisations as an ersatz of the double 
logarithmic terms ln 2 (x) resummation 

According to the philosophy of DGLAP approach, structure functions of the nucleon are a convolution of 
the coefficient functions and the evolved parton distributions. In this formalism, the polarised structure 
function gi(x, Q 2 ) for the proton is given by |13| 

g p A*,Q 2 ) = l(e 2 )[c NS ®Aq NS (x,Q 2 ) 

+ Cs®Aq s (x,Q 2 ) + 2N f C G ®AG(x,Q 2 )^, (10) 

where 

3 i=l 

Here, ei denotes the electric charge of the i quark-flavour, Aq NS , Aq s , AG are respectively the nonsinglet 
and singlet quark and the gluon polarised densities (helicity distributions). The coefficient functions d 
are computed to a given order in a 3 - PQCD evolution equations for the quarks and gluons distribution 
functions need the nonperturbative input quantities at some initial scale Q%. These input parametrisations, 
fitted to the experimental data, together with the suitable PQCD framework provide a satisfactory agreement 
of theory and measurements. The standard theoretical investigation of deep-inelastic scattering structure 
functions based on the DGLAP approach concerned originally the region of large x and large Q 2 . In this 
way the parton evolution with respect to Q 2 is taken into account, whereas the evolution with respect to 
Bjorken x is neglected. In the small- a; region, logarithms of x become also important and therefore must 
be taken into account. Assuming singular as x — > initial parton distributions ~ x~ ai (ai > 0), one can 
obtain within standard DGLAP approach a substitute of the double logarithmic In 2 (x) resummation, which 
is essential at low i < 1. However, if we take into account the double logarithmic terms via the suitable 
kernel of the evolution equations, we do not need to use the 'artificial support' in a form of the singular 
initial parametrisations. An impact of the input parton distributions on the final (after evolution) results 
is large. This is shown in Fig. 3, where we plot the LO DGLAP evolution from Qo = I GeV 2 to Q 2 = fO 
GeV 2 of the nonsinglet polarised structure function Aq NS . We test different input parametrisations of the 
general form: 

Aq NS {x,Ql)~x- a Hl-x) a \ (12) 

There is no doubt, that the small- :r behaviour of the parton densities is dominated just by the x~ ai term, 
which survives the QCD evolution when a\ > 0. Hence appropriate choice of the initial conditions must 
be consistent with used theoretical treatment. Thus there are two possible scenarios. Either we consider 
the unified evolution equations with two parts of the kernel: the standard DGLAP one and the other one 
- generating ln 2 (x) terms. Then the input distributions are assumed to be nonsingular as x — » 0. In 
this case the small- a; behaviour of the structure functions is totally governed by the evolution. Or we use 
the pure DGLAP analysis together with the singular parametrisations, which mimic the missing at low-a; 
resummation of the leading logarithms. In Fig. 4 we plot the logarithm of the polarised nonsinglet and 
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Figure 3: The LO DGLAP evolution from Q 2 , = 1 GeV 2 to Q 2 = 10 GeV 2 of the nonsinglet polarised structure 
function Aq NS as a function of x. Different input parametrisations (fT2|) . Solid: a\ — 0, a2 — 3; dashed: 
ai = 0.2, a2 = 3; dashed-dotted: a\ = 0, <i2 = 1; dotted: a\ — 0.5, a,2 = 3. 



gluon distributions evaluated at Q 2 — 10 GeV 2 within unified DGLAP+Zn 2 (a;) approach as a function of 
ln(l/x). We can estimate the effective slopes of the presented curves A(:r, Q 2 ), defined as: 



A p (,,Q 2 )= 91n [A P (*,Q 2 ) 



<91n(±' 



(13) 



Here, Ap denotes again respectively the nonsinglet quark and gluon helicity distributions (Aq NS , AG). 
From the plots we find, namely, Ajvs ~ 0.2 and Xq ~ 0.6. Hence the 'ersatz' input ~ a; _A , which is able to 
reproduce the double logarithmic ln 2 (x) resummation in the small- a; region x £ (10~ 4 ; 10" 2 ) should have 
a form: 

Aq NS ~ x-°- 2 (14) 

for the nonsinglet part and 

AG ~ x~ 6 (15) 
for the gluons respectively. As one can see, the behaviours (|14D . (|15|) are less steep than their asymptotic 
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Figure 4: The logarithm of the polarised quark nonsinglet (solid) and gluon (dotted) distributions evaluated at 
Q 2 = 10 GeV 2 within unified DGLAP+/n 2 (x) approach as a function of ln(l/x). An illustration of the slope 
A, defined in IT31). 



limits as x — > 0: 



and 



Aq NS (x -> 0) ~ g?*(x -> 0) ~ x 



NS/ 



AG(x -> 0) ~ gf(x -»• 0) ~ a: -0 ' 8 



(16) 
(17) 
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The results (|16[1 . (|17|l were obtained in [15^ via estimation of the anomalous dimensions and also in [17] - 
within IREE (infrared evolution equation) formalism. 

In conclusion, the power-like behaviour x~ x of the quark and gluon polarised distribution functions, 
generated by the resummation of the double logarithmic terms ln 2 (x), can be also obtained via the singular 
factors in the initial parton distributions. Finally, let us shortly discuss the possible values of ai in the 
leading term of the initial parton densities. The choice of the value of a\ in the input parametrisations (I12|) . 
which controls the singular small- a; behaviour, depends on the evolution length (Q 2 — Qo). If one assumes 
a very low input scale Qo < 1 GeV 2 , then already the smaller value of ai = 0.2 in the nonsinglet case 
can 'mimic' the ln 2 (x) effects. In contrast, for longer Q 2 ~ 4 GeV 2 , what denotes the shorter evolution, 
one should use more singular input with ai ~ 0.4 for the nonsinglet case. Similar (or even more singular) 
input parametrisations of the spin-dependent parton distributions have been assumed e.g. in [5], [BJ, |32) . 
In the next section we shall compare our theoretical predictions based on either the unified DGLAP+Zn 2 (:r) 
approach or the DGLAP analysis alone together with the singular inputs, with experimental data. 



5 Comparison with experimental data 

In this section we shall present our results obtained using the unified DGLAP+Zn 2 (:r) approach with 
'very running' a 3 (0. We shall also apply an alternative scenario, described in the previous section, in which 
the pure DGLAP analysis is accompanied by the singular input parton densities at the low scale Qq — 1 
GeV 2 . In this latter approach we shall compute i.a. the truncated Mellin moments of structure functions 
using directly the evolution equations for truncated moments, derived recently in [24] . Let us recall now 
some basic formulas concerning this approach. 

The evolution equations for the truncated moments of the parton densities have the form: 

dq n (x ,Q 2 ) a a (Q 2 ) , 2 

dings = -&r (p ®^ X °>Q )• ( 18 ) 

P'(n,z) = z n P(z), (19) 

where P(z) is the well-known splitting function from the DGLAP equation. q n (xo,Q 2 ) denotes the nth 
Mellin moment of the distribution function q(x,Q 2 ) truncated at xq: 



q n (x Q ,Q z ) = j dxx n ~ L q(x,Q z ). (20) 

XQ 

This formula is obviously valid also in the spin-dependent case where one replaces q by Aq, q n by Aq„ and 
P by AP e.i. the unpolarised quantities by their 'polarised' versions. It is particularly interesting to note 
that the evolution equation for the nth truncated moment has the same form as that for the parton density 
function itself with the modified splitting function P' (|19p . The truncated moments approach refers directly 
to the physical values - moments (rather than to the parton distributions), what enables one to use a wide 
range of deep-inelastic scattering data in terms of smaller number of parameters. In this way, no assumptions 
on the shape of parton distributions are needed. Using the evolution equations for the truncated moments 
one can also avoid uncertainties from the unmeasurable very small x — > region. Eq. (|18p does not account 
for the double logarithmic ln 2 (x) terms resummation, what can be mimiced by the appropriate input, as it 
was described in the previous section. This is the motivation that we use the equations for the truncated 
moments in the studies presented here. 

In order to compare theoretical predictions with experimental data over the kinematic range explored 
one should generalize the results to the small-Q 2 < 1 GeV 2 region. Thus, one can use the prescription 
introduced in [25] and applied in the studies [141 1151 [26l 127] , valid for arbitrary Q 2 : 



Q 2 ^Q 2 + Ql (21) 

and 



x = (Q 2 + Q 2 o)/(2pq). (22) 



After this rearrangement the structure function gi can be extrapolated to the low-Q 2 region (for fixed 2pq) 
including the point Q 2 = 0, although perturbative Q 2 -power and higher twist corrections may also play a 
role in this region [25]. Taking into account the small-Q 2 corrections is particularly important when one 
studies recent COMPASS measurements obtained for very small 4 • 10 -5 < x < 2.5 ■ 10 -2 at simultaneously 
very low Q 2 < 1 GeV 2 [29]. We use in our analysis input parametrisations of polarised parton densities at 
the initial scale Qq = 1 GeV 2 in a simple general form: 

Aq(x,Q 2 )=r,x- a i(l-x) a \ (23) 

where n is a normalization factor. The exponent —a\ controls the behaviour of Aq in the small- x region and 
the factor (1 — x) a2 ensures the vanishing of the parton density as x — > 1. The singular part x~ ai , where 
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Figure 5: The spin dependent structure function for proton g\ vs x, compared with HERMES data. Q 2 is the 
measured mean value (Q 2 ) at each x. Plots for different input parametrisations of the valence quarks ~ (1 — x) a2 
from up to bottom at x = 0.01: a-i = 3, a% — 2, 02 = 1. Solid (dotted) line corresponds to the positive (negative) 
solutions for polarised gluons AG . Error bars represent the statistical uncertainties. 





Q 2 


x x 


X2 


X2 

J dx gi 


Experiment 


N 


10 


0.004 


0.7 


0.04904 1 
0.04871 2 


0.051 C 
±0.003 (stat.) 
± 0.005 (syst.) 


NS 


5 


0.021 


0.9 


0.1766 a 
0.1718 b 
0.1486 c 


0.1479 H 
±0.0055 (stat.) 
± 0.0142 (syst.) 
±0.0055 (par.) 
± 0.0049 (evol.) 



Table 1: Comparison of integrals of gf = {g{ + g?)/2 and gf s with COMPASS (C) and HERMES (H) data, 
gf results for both gluon scenarios: 1 AG < and 2 AG > are shown. For g^ s the result incorporating In 2 
resummation a is compared with LO DGLAP solutions for singular input parametrisations ~ x~ ai (1 — a;) 3 with 
a x : 6 0.1-=-0.2, c 0.4. 



a\ > 0, can mimic the resummation of the leading logarithms ln 2 x. 

Figs. 5-8 and Table I contain our numerical results. Fig. 5 shows the spin dependent structure function 
for proton 

9l{x,Ql) = l(e 2 )[Aq s (x,Q 2 ) + Aq NS (x,Q 2 )] (24) 

as a function of x, compared with HERMES data [30]. Here, (e 2 ) is given by We obtain our results 

solving the unified evolution equations 0,(0 with 'flat' parametrisations ~ (1 — x)™ 2 of the parton den- 
sities. We present plots for different values of the parameter a^: 3, 2, 1 and for a negative and positive 
parametrisation of gluons. In Fig. 6 we plot 

9i =Ugf + 9i)--^h- (25) 

* 1 — jUd 

as a function of x together with COMPASS data [ST]. Here, g^, g" and gf denotes the polarised structure 
function of proton, neutron and deuteron, respectively and uid ~ 0.05 is the D-state admixture to the 
deuteron wave function. Results are shown for different contributions of gluons to the proton's spin at initial 
scale Ql, namely AG(Qo) = -0-25, 0, 0.25, 0.5, 0.75, where 



1 

AG(Q 2 ) = j dx AG(x,Q 2 ). (26) 


In Table I we collect the integrals of and s over the range of x from COMPASS and HERMES 
experiments. 
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Figure 6: = {g\ + .g™)/2 vs x, compared with COMPASS data. Q 2 is the measured mean value (Q 2 ) at each 
x. Results shown for five different contributions of gluons AG to the proton's spin: (from up to bottom) -0.25, 
0, 0.25, 0.5, 0.75. Error bars represent the statistical uncertainties. 




Figure 7: The polarised structure function g\ for the proton vs x at Q 2 = 10 GeV 2 . Results shown for five 
different contributions of gluons AG to the proton's spin: (from up to bottom) -0.25, 0, 0.25, 0.5, 0.75. 



Our results for the function as well as its first moment are in a very good agreement with the 
experimental COMPASS data. There is certain discrepancy between our predictions and HERMES data. 
This is particularly visible for the contribution to the Bjorken sum rule 



Ibjs(xi,X2, Q 



ax gi (x, Q 



dx[gl{x,Q 2 ) - g"(x,Q 2 



(27) 



shown in TABLE I. Some ansatz (input parametrisation) must be adopted in two degrees of evolution. 
Values of gi measurements in the two or three Q 2 bins for each x must be evolved to their mean Q 2 and 
then averaged. Also, the evaluation of the first moment of the structure function <?i requires the evolution 
of all measurements to a common Q 2 . In HERMES analysis this is done by using a fitted parametrisation 
[5], which increases as x -> 0: gf s ~ a;" 8 . COMPASS Group have used several fits [3 [32] which 
have been averaged. The discrepancy between our results and HERMES data reflects the fact that the 
fit used by HERMES collaboration is significantly different from ours (|23jl with a\ = 0.0. Note also that 
very close to the HERMES value for Ibjs is our result obtained within LO DGLAP approach with use of 



the singular input parametrisation gi (Qo) 



This makes the contribution from the small- x region 



< x < 0.021 more significant - at level of 30% of the total BJS compared to our estimation based on 
the unified DGLAP+Zn 2 (:r) theoretical analysis, which gives about 17%. Furthermore, it can be seen from 
Table I that LO DGLAP evolution with the appropriate input ~ x~ ai (ai > 0) for a given region of x can 
reproduce the result of the DGLAP+Zn 2 (a;) approach. In this way, a suitably chosen initial parton density 
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Figure 8: Integral of the spin dependent nonsinglet structure function g^ s = g\—g™ over the range 10 5 < x < 1 
as a function of the low- a; limit of integration. Q 2 — 10 GeV 2 . The comparison for different a\ in the input 
parametrisation ~ x~ ai (1 - x) 3 at Ql = 1 GeV 2 . Plots (from up to bottom): ai = 0, 0.2, 0.4, 0.8. 



can compensate missing low- a; effects in QCD analysis. 

We would like also to pay special attention to the evolution equation for truncated moments of the 
parton distributions. Fig. 8 illustrates the truncated contribution to the Bjorken sum as a function of the 
truncation point x. Solving the equation for moments (|18p . (|19p we test different input parametrisations 
and find the small- a; contribution 7sjs(0, 0.01, 10) (|27|) being between about 6% for the flat input ~ (1 — x) 3 
and about 60% for the very steep one ~ x~°' 8 (l — x) s . The problem of the low- a; part of the Bjorken sum 
we have also discussed in [231 133] . 

Finally, let us discuss the dependence of the polarised nucleon structure functions on the gluon distri- 
bution Ag. From Figs. 5-6 and Table I one can see that the predictions for gf and gi ~ gf (|25[) in the 
available experimentally a>region (x > 0.003) are compatible with the data independently of the assumed 
gluon function. Large experimental uncertainties for low x do not allow one to discriminate between differ- 
ent, in particular positive and negative polarised gluon densities. In Fig. 7 we compare the proton structure 
function g\ at Q 2 = 10 GeV for different fractions of the nucleon spin carried by gluons at the initial scale 
Qq. Note, that g\ essentially depends on the gluon distribution only for very low x - not before x « 0.01. 
Our parametrisations of AG (|26p reflect the latest experimental determinations of the gluon polarisation 
at COMPASS [32], RHIC [35] and STAR [3B]. The shape of AG(x,Q 2 ) is poorly known and the present 
experimental data support both positive and negative distributions, resulting in small | AG |~ 0.2 to 0.3 
(COMPASS) or large AG = -0.56T2.16 (RHIC), AG = -0.45 to 0.7 (STAR). It is possible that a significant 
contribution to AG comes from low- a;. A knowledge of the small- a: behaviour of AG(a;,Q 2 ) would provide a 
constraint on the shape and the sign of the gluon component. We hope future measurements at RHIC over 
a wide range of x and Q 2 will enable precise determination of the gluon contribution to the nucleon spin. 

6 Conclusions 

In this paper we have presented results for the spin structure function g\ of the nucleon together with 
comparison with latest HERMES and COMPASS data. We have applied an approach that combines LO 
DGLAP Q 2 evolution with the resummation of the double logarithmic terms ln 2 (x). This unified framework 
goes beyond the standard LO and NLO PQCD evolution of the parton densities and becomes essential for 
x < 0.01. In our analysis, we have focused on the taking into account the 'very running' coupling effects. For 
the small- a; region the more justified is the use of a s — a s (Q 2 jz] instead of a s — a s (Q 2 ), with z being the 
longitudinal momentum fraction of a parent parton, carried by a next generation. In this way a s becomes 
'very running' i.e. runs in each ladder rung depending on the gluon virtuality. We have shown, that the 
impact of these running coupling effects becomes important at x < 0.01 and significantly damp the results. 
The decreasing factor at very small x = 10 -5 can be about 1/2 for the nonsinglet and 1/6 for the singlet 
(gluon) distribution function in comparison to the standard a s (Q 2 ) prescription. 

In order to calculate the first moment of gi over the available experimentally x region, we have solved 
the direct evolution equations for truncated moments of the parton densities. In this approach we have 
utilized the fact that the resummation of the double logarithmic terms, missing in the standard DGLAP 
approximation, can be mimiced by singular input parametrisation of the parton distributions. The truncated 
moments approach refers to the physical values - moments (rather than to the parton distributions), what 
in future analyses could enable one to use a wide range of deep-inelastic scattering data in terms of smaller 
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number of parameters. In this way, no assumptions on the shape of parton distributions are needed. 

Our theoretical predictions for the polarised structure function <?i and its first moment for the deuteron 
are in a very good agreement with COMPASS data. There is certain discrepancy between our predictions 
and HERMES data, particularly visible for the contribution to the Bjorken sum. This reflects the fact that 
the fit used by HERMES collaboration is significantly different from ours. It must be emphasized, that the 
final (after evolution) results strongly depend on the the input parton distributions assumed. 

Finally, let us discuss the dependence of the polarised nucleon structure functions on the gluon distri- 
bution Ag. Large experimental uncertainties in the low- a; region do not allow one to discriminate between 
different polarised gluon densities. The shape of AG(x,Q 2 ) is poorly known and the present experimental 
data support both positive and negative gluon distributions. A knowledge of the small- £ behaviour of the 
gluon component and possibly a significant contribution to AG from this region would enable to resolve the 
nucleon spin puzzle. This is a challenge for future theoretical and experimental efforts. 
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